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Abstract 

In this paper, the decay constants and mean square radii of pseudoscalar heavy mesons are 
studied in the SU{3) symmetry breaking. Within the light-front framework, the ratios /d s //d and 
fB a /fB are individually estimated using the hyperfme splittings in the D* s ^ — and B*, \ — B^ 
states and the light quark masses, m s ^ q {q = u,d), to extract the wave function parameter /3. 
The values /d 3 //d = 1-29 ± 0.07 and fsj fs = 1-32 ± 0.08 are obtained, which are not only 
chiefly determined by the ratio of light quark masses m s /m q , but also insensitive to the heavy 
quark masses m c ^ and the decay constants fr>,B- The dependence of fsd 1 fs on AMb c b* with 
the varied charm quark masses is also shown. In addition, the mean square radii are estimated 



as well. The values ^ '(r£+> / '(r 2 D+ ) = 0.740+^ and ^ (r£j> / (r% ) = 0.711^;^ are obtained, 
and the sensitivities of (rp) on the heavy and light quark masses are similar to those of the decay 
constants. 
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I. INTRODUCTION 



The decay constants of pseudoscalar heavy mesons with c and b quarks play an important 
role for studies of CP violation and in extracting the Cabibbo-Kobayashi-Maskawa matrix 
elements. Experimentally, new data on the charm meson decay constants fjj and fo s have 
been reported [l, 2\. As the calculations of the decay constants are related to the wave func- 
tion overlap of the quark and antiquark which are governed by the strong interaction, they 
therefore provide a crucial manner to compare different theoretical methods. In addition, 
the determination of fs s remains beyond the reach of current experiments, thus the reliabil- 
ity of estimated fs s by a theoretical approach is dependent on whether the determinations 
of fr> and fr> g by this approach are consistent with the new data. During the last decade, 
;he decay constants of pseudoscalar heavy mesons have been studied in lattice simulations 



9|, in the QCD sum rules approach |lOHl3j. and in the relativistic quark model |l4j-|l8|. 
The understanding of the electromagnetic (EM) properties of hadrons is also an impor- 
tant topic, and the EM form factors which are calculated using nonperturbative methods 



are the useful tool for this purpose. There have been numerous experimental 
theoretical studies 



19 



-241 and 



25] of the EM form factors of the light pseudoscalar meson (ir and K). 



However, the 
fewer studies 



M form factors of heavy mesons (which contain one heavy quark) have much 
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271 ] than those of light ones. The present paper is devoted to an analysis 



of the wave function and decay constant by the hyperfme mass splitting of heavy mesons 
and the formulas of the decay constant and mean square radius within the light-front (LF) 
framework. We present the SU(3) symmetry breaking effect in decay constants and electro- 
magnetic properties of pseudoscalar heavy mesons. 

The light-front quark model (LFQM) is a promising analytic method for solving the 



nonperturbative problems of hadron physics 



28 



as well as offering many insights into the 



internal structures of bound states. The basic ingredient in LFQM is the relativistic hadron 
wave function which generalizes distribution amplitudes by including transverse momen- 
tum distributions and contains all the information of a hadron from its constituents. The 
hadronic quantities are represented by the overlap of wave functions and can be derived 
in principle. The light-front wave function is manifestly a Lorentz invariant, expressed in 
terms of internal momentum fraction variables which are independent of the total hadron 
momentum. Moreover, the fully relativistic treatment of quark spins and center-of-mass 



2 



motion can be carried out using the so-called Melosh rotation 29] . This treatment has been 
successfully applied to calculate phenomenologically many important meson decay constants 



and hadronic form factors 
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36|. 



The remainder of this paper is organized as follows. In Sec. II an analysis of wave 
function and decay constant is presented. In Sec. Ill the formulism of LFQM is reviewed 
briefly, and the formulae of decay constant and mean square radius are derived. In Sec. IV 
numerical results and discussions are presented. Finally, the conclusions are given in Sec. 
V. 



II. ANALYSES OF WAVE FUNCTION AND DECAY CONSTANT 

The decay constant fp for a pseudoscalar meson is defined by a matrix element of the 
axial vector current between the vacuum and the meson bound state: 

(0\qi W q2\P(P))=ifpP„. (2.1) 
In a nonrelativistic approximation, fp is related to the Bethe-Salpeter wave function at the 



origin |\P(0)| as 



37 



39] 



h - $gl*(0)|, (2.2) 

where N c is the color number and Mp is the mass of the meson. If we consider the potential 
of a hyperfine interaction inside the meson to 0(a s ): 



Vhf = ^ § + — o 8*{r), (2.3) 



where a s is the strong coupling constant, £1,2(^1,2) are the spins (masses) of the constituent 
quark. For the s-wave meson, the first term of Eq. (12.31) has no contribution and the second 
term can distinguish the pseudoscalar and vector mesons. Therefore, the hyperfine mass 
splitting is obtained as 

AM Py = |^|*(0)| 2 . (2.4) 
9mim2 

By combining Eqs. (12. 2p and (I2.4p and canceling |\l/(0)|, we obtain: 

f27AM PV m im2 \ 1/2 
fp ~ { 8na s M P ) ' (2 ' 5) 



If we suppose the strong coupling constants a s (D) ~ a s (D s ) and a s (B) ~ a s (B s ), then the 
ratios of the decay constants can be obtained as 

/d. = ( AM DsDt M D m s \ 1/2 /b. = / AM Bs b* M b m s \ 1/2 

where q = u,d. From Eq. (12.6p . we find that the ratios are dependent on the ratio of 
light quark masses m s /m q and are independent of heavy quark masses. Furthermore, by 
canceling the ratio m s /m q , we have a relation which does not contain any parameter in the 
nonrelativistic approximation: 



f Bs (AM DD .AM BaB . M Ds Mb\ 1 ' 2 f Da 



(2.7) 



f B \AM DsD *AM BB * M d MbJ fn 

If one wants to include the relativistic correction to the ratios of the decay constants, 
not only the values of m Siq , but also the form of wave function ty(f) must be known. Let us 
come back to Eq. ( 12. 2(1 . The deviation of Eq. ( 12.2(1 uses the Fourier transform 

* (f > = / J0r^ k ^ ^ 

where <f>{k) is the wave function in the momentum space. If the Fourier transform is evaluated 
for the positive-energy projection of the Bethe-Salpeter wave function at equal "time" z + = 



z° + z 3 = 0, then 



40] 



where x is the longitudinal momentum fraction, k± are the relative transverse momenta, 
and <j)(x, k±) satisfies the normalization condition: 

dxd 2 k i 



{x,k ± )\ 2 = l. (2.10) 



2(2vr) 3 

In general, the momentum distribution amplitude <p(x, k±) is obtained by solving the light- 
front QCD bound state equation Hlf\P) = M\P) which is the familiar Schrodinger equation 
in ordinary quantum mechanics, and Hlf is the light-front Hamiltonian. However, at the 
present time, how to solve the bound state equation is still unknown. Alternatively, we 
come back to the still- unknown wave function \I/(f) and express it as a linear combination 
of the arbitrary known functions which form a complete set. Of course, the complete set is 



4 



not unique. Here we give two examples for comparison. One is the solution of 1/r potential 



4l|: 



Kim(r) 



. n 



3/2 



(n-l-l)\ 



Ln-U^r/n)e W (-^) Y lm , (2.11) 



_2n[(n + /)!] 3 

where superscript c means "Coulomb," (3 is a parameter which has the energy dimension, 
Yi m is the spherical harmonics, and L q _ p (x) is an associated Laguerre polynomial which is 
defined as 



T p 

^q-p 



e~ x x q ) 



(2.12) 



The other complete set is the solution of an isotropic harmonic oscillator 411 ]: 



/3 3 / 2 , / B 2 r 2 \ 

where superscript g means Gaussian and the first few h n i(j3rys are 

~ 4 



(2.13) 



h 



oo 



2, ^ii 



3' k22 
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2() 



>/6 1 - ~/?V 



(2.14) 



Then the exact solution can be expressed as 



^(r) = E4M(r), 



(2.15) 



where X] 



c (g) i2 



nim \ u 'nlm 



1. This way seems very clumsy because, apart from /3, it also in- 



troduces a series of undetermined coefficients, a n i m . The following considerations, however, 
improve the situation. First, only the coefficients a n00 survive because we just studied the 
s-wave meson. Second, we substitute Eq. ( 12. 15}) to Eq. ( 12. 41) and obtain 



AM, 



327T«, 



0" 



PV 



/] a nOO^no(0) 



(2.16) 



9mim2 \Att 3 / 2 t 

which takes the Gaussian case, for example. It is worth noting that the square bracket in 
Eq. (12.161) is independent of (3. Then, the ratio of hyperfine mass splittings can be reduced 
as: 

AM DlD . m„ f (3 CS \ 3 AM BsBt 



m q ( f3 cs 

m s \(3cq, 



m q f (3 bs 



(2.17) 



AM Dfl . m s \/3 cq J AM BB . 

which does not include any coefficient a n j m . Equation (I2.17P is also suitable to the Coulomb 
case. In fact, due to the 5 function in Eq. (12. 3p having the dimension of an energy cube and 



r*is vanishing here, Eq. (12. IT)) holds for any wave function which contains only one hadronic 
parameter (3. In addition, the ratios (3 CS / (3 cq and f3bs/Pbq in Eq- (12.1 T)) are mainly influenced 
by the ratio m s /m q . This situation leads to the ratios /dJId an d /b s //b are sensitive to 
the SU(3) symmetry breaking, but are insensitive to the heavy quark masses, which will be 



shown later. In the literature, there are some early attempts [42|, |43( to account for flavor 
symmetry breaking in pseudoscalar meson decay constants. In the next section the wave 
function $?(r) and the values of m SA are studied within the light-front framework. 

III. LIGHT-FRONT FRAMEWORK 
A. General Formulism 

An s-wave meson bound state, consisting of a quark qi and an antiquark q 2 with total 



momentum P and spin J, can be written as (see, for example 32]) 



\M(P, S, S z )) = J {d 3 Pl }{d 3 P2 } 2(2tt) 3 5 3 (P -px - p 2 ) 



x $ ss *(pi,p 2 ,\ 1 ,\ 2 ) |gi(pi,Ai)g 2 (p 2 ,A 2 )), (3.1) 

Ai,A2 



where p\ and p 2 are the on-mass-shell light-front momenta, 



V = (p + ,P±) , P± = (p\p 2 ) , P~ = q ■ , (3-2) 



P 



and 



{d3p} ^ dp+fpj. 



2(2tt) 3 ' 

|g(pi,A 1 )g(p 2 ,A 2 )) = 6t(p 1 ,Ai)dt(p 2 ,A 2 )|0), (3.3) 
W, A'), b\p, A)} = {d(p r , A'), d\p, A)} = 2(2tt) 3 5 3 (p' - p) 5 yx . 

In terms of the light-front relative momentum variables (x, k±) defined by 

p+ = (l-x)P + , p+ = xP + , 

Pi± = (1 - x)P ± + k ± , p 21 _ = xP ± - k ± , (3.4) 
the momentum-space wave-function fy SSz can be expressed as 

® ss *(pi,p 2 ,\i,\2) = -^fR S x%(u,K ± ) <P(x,k ± ), (3.5) 

V Mr. 



where Rf^L constructs a state of definite spin (S,S Z ) out of light-front helicity (Ai,A 2 ) 
eigenstates. Explicitly, 

11 11 

R \?Za( x > k ±) = J2 ( A iI 7 ^m( 1 ~ x,k±,m 1 )\s 1 )(\ 2 \'Jl\ 1 (x, -k±,m 2 )\s 2 )(--;s 1 s 2 \--;SSj^.6) 



where |sj) are the usual Pauli spinors, and TZm is the Melosh transformation operator 30l.l3lj: 

rrii + XiM + ia s \ ■ k± x n 



(s\n M (xi,k±,mi)\X) 



(rm + x l M ) 2 + k\ 

with x\ = 1 — x, x 2 = x, and n = (0, 0, 1) as a unit vector in the z direction. In addition 



(3.7) 



M, 



ei + e 2 2 = -4 - + ' 



1 — x 



X 



(3.1 



+ k 2 , + kl 



m" 



where k z is the relative momentum in z direction and can be written as 



xM m\ + k\ 



(3.9) 



2 2xM 

Mq is the invariant mass of qq and generally different from the mass M of the meson which 
satisfies M 2 = P 2 . This is due to the fact that the meson, quark and antiquark cannot be 
simultaneously onshell. We normalized the meson state as 

(M(P', S', S' Z )\M(P, S, S z )) = 2(2n) 3 P + 5 3 (P> - F) Ws'S, , (3.10) 



which led to Eq. ( |2~T0|) . 



ss z 



In practice, it is more convenient to use the covariant form of F Ai " 



30 



31 
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y/2 M (M + m l + m 2 



-u(pi, A] 



M )Tv(p 2 ,X 2 ), 



ii- 
(3.11) 



where 



M = ^M 2 - (mi - m 2 ) 2 , 

u(p, X)u(p, A') = ^5 X ,\>, 
p+ 

v(p,X)v(p, X') = -^S x> x', 
p + 



P=Pl +P2, 

J2u(p,X)u(p,X) = 

A 

J2 V (P, x )v(p, A) 

A 



p 1 + m 



P 



m 



P 



(3.12) 



For the pseudoscalar meson, we have T = 75, Eq. f)3.1ip can then be further reduced by the 
applications of equations of motion on spinors 



Pi P2 
V2M 

7 



u{pi,Xl)l5V(p2,X 2 ). 



(3.13) 



Next, we derive the formulas of the decay constant and the mean square radius for the 
pseudoscalar meson. The former is the main subject of this work, and the latter is used to 
fix some parameters. 



B. Formulas for decay constant and mean square radius 

The decay constants of pseudoscalar mesons P(qiq~2) are defined in Eq. (12. ip . The matrix 
element can be calculated using the formulism in the last subsection: 

m2l»lMi\P{P)) = /{rf 3 Pi}{^ 2 }2(27r) 3 5 3 (P-p 1 -p 2 )0 P (x,fc ± )^ ( ; A2 (x,fc ± ) 

x (0|<? 2 7M75gi|M2>. (3.14) 



Since M yx(l — x) = J A 2 + kj_, the decay constant can be extracted as: 

.4 



f P = 2^2N C f {dx} P (x, k ± ) (3.15) 

J ■ * A 2 + k 2 



where {dx} = dx ^^ and A = vti\x + m 2 (l - x). 

Next, the EM form factor of a meson P is determined by the scattering of one virtual 
photon and one meson. It describes the deviation from the pointlike structure of the meson, 
and is a function of Q 2 . Here, we considered the momentum of the virtual photon in a 
spacelike region, so it was always possible to orient the axes in such a manner that Q + = 
(P' — P) + = 0. Thus, the EM form factor was determined by the matrix element: 

(P(P')\J + \P(P)) = e F P (Q 2 )(P + P') + , (3.16) 

where J M = qe q e^q, e q is the charge of quark q in e unit, and Q 2 = (P' — P) 2 < 0. With 
the light-front framework, Fp can be extracted by Eq. ( 13.161) : 

/A 2 -\- k± ■ k' 
{dx} - 1 =<j) P (x, fc_i_)0p'(x, k' ± ) 

y/ A 2 + k 2 ± y/ A 2 + kf 

+ e q - 2 {dx}^= — -, X =(/)p(x,k ± )(l)p'(x,k'l), (3.17) 
I A 2 + k 2 J A 2 + k'[ 2 



where k' ± = k± + xQ±, k" L = k± — (l — x)Q±. For applying this to Eq. (13. 22 p . it is convenient 



to consider the term 0p = 4>p(x, k±)/ \J A 2 + k\ and take the Taylor expansion around k\ 

a2 - d 2 



4>p,(k'l) = 0p,(kl) + 



dki 



(k't-kiY + (3.18) 

o 



Then, by using the identity 

J d 2 k ± {k x ■ A ± )(k± -B L ) = ^J d 2 k ± k\ A ± ■ B A 
we can rewrite (13.171) to 



(3.19) 



Fp{Q 2 ) 



[e qi + e q2j 

l2 [ r. i ±1 



— Q J {dx} 4> P (x, k±) 
+ 0(Q 4 ), 



x 2 e qi + (1 - x) 2 e 



in 



e, 



where 



d 2 4> M " 

,2 \2 



(3.20) 



(3.21) 



4>m V dki J 1 4>m\ 

It should be realized that the size and the density of a hadron depend on the probe. For an 
EM probe, it is the electric charge radius (r 2 ) 1 / 2 that is obtained. In the experimental view, 
(rf>) cannot be measured directly and is obtained by fitting the slope of Fp(Q 2 ) at Q 2 = 0, 



i.e., 



(r%) = 6 



dF P (Q 2 



dQ* 



(3.22) 



Here the mean square radius is easily obtained: 



(r P ) = K) + (r|) 

= e 91 | — 6 J {dx}x 2 (f)p 



(A 2 + 2kl)4r + (A 2 + kl)k] 



2i 



dk 



dki 



+ e, 



<-n 



6 {dx}(l- xf(\) 1 



{A 2 + 2kl)4r + {A 2 + kl)k 2 ^ 1 



2-\ 



dk\ 



dk 2 ±/ 



•(3.23) 



It is worth mentioning that, first, the static property -Fp(O) = ep is quite easily checked in 
Eq. (13.201) . Second, from Eq. (13.231) . we find that the mean square radius is related to the 
first and second longitudinal momentum square derivatives of (f> which contain the Melosh 
transformation effect. 

If we take the heavy quark limit ni\ = rriQ — > oo, rriQ(Mp) is unimportant for the low 
energy properties of the meson state, so it is more natural to use velocity v instead of 
momentum variable P. The normalization of the meson state is rewritten as 45] 



(P(v')\P(v)) = 2(2it) 3 v + 5 3 (Av - At/), 



(3.24) 



where A = Mp — tuq is the residual center mass of the heavy meson and the meson states 
have a relation \P(v)) = (Mp)~ 1 ^ 2 \P(P)) . In addition, since x is the longitudinal momentum 
fraction carried by the light antiquark, the meson wave function should be sharply peaked 
near x ~ Agc^/mg. It is thus clear that x — >■ and only terms of the form X = xttiq 
survive in the wave function as mq — > oo; that is, X is independent of mq in the heavy 
quark limit. Therefore, the normalization of the wave function Eq. (I2.10p is rewritten as 

dXd 2 k\ 



MX,k±)\ 2 



1. 



(3.25) 



where <f(X, kj 



l m Qj 



2(2vr) 3 

x,kj_). Other replacements are A — >• A = X + m q2 and 



(x, k±) —> (p(X, k±). Thus we can rewrite (13.231) as 



(r Qq2 ) = (r 2 Q ) + (r|) 

= -^-i- 6 f{dX}X' 
m Q «■ J 



(A 2 + 2kl)^L + (A 2 + kl)kli 



+ e. 



92 



6 / {dX}(p 



(A 2 + 2k] 



^ + (A 2 + k 2 ± )kl\ 



d 



dk\ 

\ 2i 



dkl t 



(3.26) 



The first term of Eq. (I3.26|) vanished when itlq — > oo. This means that not only (^q 92 ) is 
blind to the flavor of Q, but also (rf>) is insensitive to mi for the heavy meson. The former 
is the so-called flavor symmetry and the latter will be proven in the numerical calculation. 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

In the nonrelativistic (NR) approximation, we substituted the experimental data [l| to 
Eq. (12. 7L and obtained /b s //b = (1.03 ± 0.02) fpj fr>. If one wanted to evaluate fsj fs 
and foslfo individually, then m s = 483 MeV and m„ = 310 MeV were the "best-fit" values 
for the pseudoscalar and vector light meson masses 46[. The ratios were 

NR fB 

The former was a little smaller than the data 

M fDjf D \e, P 

was almost larger than the other theoretical calculations (see Tables II and III). 

In the light-front framework, the momentum distribution amplitude <fi(x, k±) or the wave 
function ty(r) in principle is unknown unless all the coefficients a n oo are obtained. However, 
we may suppose a° 00 = 1 or clq 00 = 1, that is, \&(f) = ^i 00 (r) or ^(r) = ^Q 00 (r) as a trial 



fDs 

fn 



1.226 ±0.002, 



NR 



1.24 ±0.02. (4.1) 
1.27 ±0.06, and the latter 



10 



wave function to fit the relevant data. The other a n oo's will be subsumed if the parameters 
appearing in the momentum distribution amplitude can not satisfy all experimental results. 
In other word, the coefficients a n i m are taken as another kind of parameter. Of course, based 
on the principle of quantum mechanics, the physical meanings of these new parameters are 



clear. Here we list the first ^rm'- 



(f> c 100 (x,k ± ) 



^000 ( x > 



2vr\ 1 / 2 

WJ 1 

7T \ 3 /4 
J 2 



eie 2 



x(l - x)M 



eie 2 



k\ + k 2 + f3 2 
k i ± k v 



exp 



2/3 2 



(4.2) 
(4.3) 



x(l - x)M 

and use the experimental data of f n + = 130.4 ± 0.2 MeV and (r 2 ^) 1 / 2 = .672 ± 0.008 
fm to fit the parameters m q and (3 qq . The results are m q = 0.172(0.251) GeV and (3 qq = 
0.555 =F 0.011(0.317 =F 0.007) GeV for 0f oo (0ooo)- A s for the strange quark mass, in Ref. 
47 1 m s — m u = 0.23 GeV was obtained with some interaction potentials, while in Ref. 
3l| m s — m u = 0.12 GeV in the invariant meson mass scheme. So here we use the values 
m s - m u = 0.180 ± 0.050 GeV and f K + = 155.5 ± 0.8 MeV to fix (3 sq . The results are 

i. The charge radius (r 2 <+ ) 1 ^ 2 and the mean 



aools^ooo 



$ sq = 0.463;°:^(0.354;°;^) GeV for 
square radius (r 2 K0 ) were calculated by these parameters and were listed in Table I. We 

TABLE I: Charge radius (r| r+ ) 1 / 2 and the mean square radius (f^o) of the experiment, this work, 
and other theoretical estimations. DS is Dyson-Schwinger equations; VMDx is vector meson 
dominance plus an effective chiral theory; BS is Bethe-Salpeter equation. 





Experiment [1] 


0100 


^000 


DS[48] 


VMDx [49] 


P QCD[50] 


BS [51] 


(r 2 ^) 1 / 2 (fm) 
<4o> (fm 2 ) 


0.560 ±0.031 
-0.077 ±0.010 


n 7ln +0.033 

_ n 121 -0.036 
u - izi +0.039 


u - ou ' -0.012 

-n n72-°- 017 

u - u ' z +0.019 


0.49 
-0.020 


0.616 
0.057 


0.570 
-0.0736 


0.62 
-0.085 



found that, on one hand, the value of (rj <+ ) 1 ^ 2 ( (r 2 KQ )) for 05oo was too large (small) than 
that obtained in the experiment. Then, the coefficients a c nQQ for n > 1 may be taken as 
nonzero to correct the fitting of (r 2 K+ ) 1 / 2 and (r 2 K0 ). However, the mean square radii of <^ 00 
is greater when n is larger, or (^ 2 )^ 00 < ( r2 )^ 00 < (^^ioo < ••• This means, for decreasing 
the value of (r|, + ) 1//2 , the values of a c n00 must be artificially arranged in order to cancel out 
the contributions of <^ 00 (n > 2) mutually. It is too hard to achieve now. On the other 
hand, the results for 0qoo were consistent with the experimental data. Therefore we only 
use the Gaussian-type wave function, 0ooo; to the following calculation. By combining Eq. 
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FIG. 1: Dependence of Jds/Id on fjy with m c = 1.2, 1.5, 1.8 GeV. The left and right vertical dash 
lines correspond to the lower and upper limits of the data fu = 205.8 ± 8.9 MeV, respectively. 



a. 



(12.171) . the experimental data [1], and the light quark mass in above, we obtained the ratios 

as 



Pes 

(3 cq 



1.20 ±0.04, 



1.20 ±0.05. 



(4.4) 



Obviously, the SU(3) symmetry breaking is the major contribution to the ratios in Eq. 

For the heavy quark masses, the quite different values were also used in the model cal- 
culations. For example, m c = 1.38 GeV and m& = 4.76 GeV which were fitted for the 



spectrum of the p-wave charmonium and bottomonium states 52j; and m c = 1.8 GeV and 
rrih — 5.2 GeV which were obtained from the potential models and the variational principle 
531 ] . Here the values m c = 1.2, 1.5, 1.8 GeV and m\> = 4.2, 4.7, 5.2 GeV were taken into 
account. By combining Eqs. (I3.15p . (14. 3p . (14. 4p . and the quark masses m q ( s ) = 0.251(0.431) 
GeV, the dependences of /d s //d on f D with three different m c 's and fsj on f B with 
three different m^'s were shown in Fig.l and Fig.2, respectively. It was easily found that the 
ratios jo s l Jd an d IbJ Ib were not only insensitive to the heavy quark masses m c and m b , 
but also insensitive to the decay constants /o and fs, respectively. 

Recently, the CLEO collaboration updated their data about the branching fraction for 
the purely leptonic decay D + — > and reported [l| = 205.8 ± 8.9 MeV. By using 
this value, we determined the ratio fnj /d = 1-29 ± 0.07 and the decay constant fo s = 
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FIG. 2: Dependence of /s s //s on with m c = 4.2, 4.7, 5.2 GeV. The left and right vertical dash 
lines correspond to the lower and upper limits of the data fg = 204 ± 31 MeV, respectively. 

264.5 ± 17.5 MeV with m c = 1.5 GeV. We found these results were consistent with the 
data [2]: = 261.2 ± 6.9 MeV and f%?/f%? = 1.27 ± 0.06, which were the average 

of the CLEO and Belle results (which included the radiative corrections). In addition, 
our values were generally larger than the other theoretical calculations. Table II compares 



the theoretical calculations with experimental value. For the bottom sector, the Belle 
and Babar 
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56] collaborations found evidence for B — )■ r v decay which was not 



helicity suppressed. However, the Belle and Babar values had 3.5 and 2.6 standard-deviation 
significances, respectively; thus the average was provisional [57]: B(B~ — >■ r~v) — (1.42 ± 
0.43) x 10~ 4 . We extracted the decay constant f B xp = 204 ± 31 MeV. By using this value, 
the ratio f B JfB = 1-32 ± 0.08 and the decay constant f Bs = 270.0 ± 42.8 MeV with 
m& = 4.7 GeV were obtained. Table III compares the theoretical calculations with the 
experimental value. Similar to the charm sector, our ratio /b„//b was almost larger than 
all other calculations. It is worth mentioning that the decay constants of both pseudoscalar 
and vector heavy mesons have already been investigated by the author of Ref. 18j with 



the analysis of magnetic dipole decays o 



quark model. The parameters in Ref. 18 



various heavy flavored mesons in the light-front 
were constrained by the variational principle for 



the QCD-motivated effective Hamiltonian. Roughly speaking, our above results showed the 
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TABLE II: Theoretical calculations of the decay constants fjj, fo s (MeV), and the ratio /d s //d- 
QL is quenched lattice calculations, BS is Bethe-Salpeter equation, Linear and HO are the different 
potentials within LFQM. We have quoted only the value with m c = 1.5 GeV in this work (LF). 





fD 


fn 


Id s 1 Id 


Exp er iment 


205 8 ± 8 


9 [11 


261.2 ± 6.9 [2] 


1 27 dz 06 a 

-i • i — i — V/ • v / \y 




This work (LF) 


205 8 ± 8 

*i VJ fJ . U _L_ L? 


9 


264.5 ± 17.5 


1 29 ± 07 




This work (NR) 








i oor _i_ n nno 

_i_.__.vj _i— vj • vjvj— 1 




Lattice fHPOCD+UKOCD) [31 


208 ± 4 




241 ± 3 


i i «o + 009 

1 • J_ V' — ■ I VJ • VJVJ tJ 




OL fOCDSF) [41 


«i VJ VJ i VJ _l_ 


3 ± 22 


220 ± 6 ± 5 db 11 


1 ORQ _i_ n ni o 

l.VJVJU _l_ VJ.VJIO 


+ 090 

_l_ VJ.VJ— -VJ 


OL (Taiwan) [51 


235 ± 8 ± 


14 


266 ± 10 ± 18 


1 13 dz 03 dz 


0.05 


Lattice (FNAL+MILC+HPQCD) [6] 


201 ± 3 ± 


1 7 

1 f 


249 ± 3 d= 16 


1 . _ 1 I VJ.VJ-J-. I 


07 

VJ . VJ 1 


OL fUKOCD") [81 


oi pi _i_ i rH 

_■ JL VJ —I— ± VJ_ 


-17 
-16 


238 ± 8 + ;I 


i i q _i_ n n9+° 

X.JLO _1_ VJ.VJ-__n 


04 
02 


OL [91 


911 -J- 

ill 14: 


-0 
-12 


231 ± 12+^ 


i i n _l n no 

1.1U VJ.VJi 




OOD Sum Rules [101 

V.. 1 J vJ Villi ILLIIVjO 1VJ 


177 ±21 




one; _i_ 99 


1.16 dz 0.01 dz 


0.02 


QCD Sum Rules [12] 


203 ± 23 




235 ± 24 


1.15 dz 0.04 




Field Correlators [14] 


210 ± 10 




260 ± 10 


1.24 dz 0.04 




Potential Model [15] 


234 




268 


1.15 




BS [16] 


230 ± 25 




248 ± 27 


1.08 dz 0.01 




BS [17] 


238 




241 


1.01 




Linear{HO} [18] 


211{194} 




248{233} 


1.18{1.20} 




"This value is obtained by combining /d = 


= 205.8 ±8 


.9 MeV 


Q] and f Ds = 261. 


2 it 6.9 MeV U 





flavor SU(3) symmetry breaking m s /m q = 1.72 ± 0.20 leads into the ratios 

^ = 1.29 ±0.07, 

JD 

fp- = 1.32 ±0.08. 

JB 

For the B c meson, however, both the decay constant and the hyperfine splitting have not 
been measured yet. We considered the ratio of hyperfine mass differences: 

AM BcB * m q fp bc V 
AM BB * m c \p b J 
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TABLE III: Theoretical calculations of the decay constants Jb, Jb s (MeV), and the ratio Jb„/ Jb- 
Only the value with = 4.7 GeV has been quoted in this work (LF). 





J B 


Jb s 


J Bsl J B 


Experiment 


204 ± 31 a 








This work (LF) 


204 ± 31 


270.0 ± 42.8 




1.32 ± 0.08 


This work fNR) 








1 24 + 02 

1 > JL 1 V J • V ) w 


OL fOCDSF) [41 


190 ± 8 ± 23 ± 25 


205 ± 7 ± 26 ± 17 


1 080 ± 028 ± 031 

1 • UUu — 1 — • ^ W — 1 — KJ • KJtJ 1 


Lattice fHPOCD) [71 


216 ±9±19±4±6 

1 V ' 1 kJ 1 1 ' ' 1 X 1 v / 


259 ± 32 




1 20 ± 03 ± 01 

1 • _ V / 1 V J • V J • J 1 V 7 • V / 1 


OL fUKOCDl [81 


177 ± 17 ± 22 

i i i i i i i _■ ._ 


204 ± 12±J^ 




1 15 ± 


OL [91 


179 ± 18 + n 6 

1 ) l/ — 1 — 1U Q 


204 ± 16lo 8 




1 14 ± 03 + 2' S? 


OCD Sum Rules [111 

Vc£ V_ ■ -1 V k> LI. 1.1 J. -L Kj HI V^O 11 


178 ± 14 


200 ± 14 




1 12 + 01 + 03 

1 ■ 1 1~ 1 V/ • V j 1 1 V J > \ J KJ 


OCD Sum Rules [121 


203 ± 23 


236 ± 30 




1 16 ± 05 

X • ±. KJ —i— KJ.KJtJ 


QCD Sum Rules [13] 


210 ±19 


244 ± 21 




1.16 


Field Correlators [14] 


182 ±8 


216 ±8 




1.19 ±0.03 


Potential Model [15] 


189 


218 




1.15 


BS [16] 


196 ± 29 


216 ± 32 




1.10 ±0.01 


BS [17] 


193 


195 




1.01 


Linear{HO} [18] 


189{180} 


234{237} 




1.24(1.32} 



a This value is extracted by the branching ratio: B(B — > t v) = (1.42 ± 0.43) x 10" 
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Similar to the above cases, the ratio JbJ Jb was insensitive to the value of m& and sensitive 
to that of m c /m q . The dependences of JbJ Sb on AMb c b* with m c = 1.2, 1.5, 1.8 GeV and 
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58- 



rrih = 4.7 GeV were shown in Fig.3. Some model predictions were made for Jb, 
621 ]. and the range of these values was J Bc = 360 ~ 517 MeV or JbJ Jb — 1-76 ~ 2.53. As 
shown in Fig.3, this range corresponded to AM BcB * = 16 ~ 64 MeV. We found that this 



result was consistent with a calculation using the nonrelativistic renormalization group 63] 
AM Bc b* = 48 ± 15t\i MeV. 

Besides, the mean square radii of the heavy meson are calculated by the above parameters 
and Eq. (13 .23j) . The dependences of (r 2 D+ D0 D ) on m c and (r B+ B0 B ) on m& were shown in 
Fig. 4 and Fig. 5, respectively. It was easily found that, as mentioned at the end of sec. Ill, 
the mean square radii (r 2 D+ D0 Ds ) and (r B+ B0 Bs ) were insensitive to the heavy quark masses 
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FIG. 3: Dependences of f B Jf B on AM BcB * with m c = 1.2, 1.5, 1.8 GeV and m b = 4.7 GeV. The 
low and high horizontal dash lines correspond to f Bc = 360 MeV and f Bc = 517 MeV, respectively. 
The left and right vertical dash lines correspond to AM BcB * = 16 MeV and AM BcB * = 64 MeV, 
respectively. 
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FIG. 4: Dependences of {r 2 D+ D o D ) on m c = 1.2 ~ 1.8 GeV. 



m c and m&, respectively. We used m c = 1.5 GeV and mb = 4.7 GeV to estimate the mean 
square radii of the heavy meson and the results are listed in Table IV. It is interesting to 
note that the value (r B+ ) is slightly lower but comparable to what one would obtain from 



the lattice calculation of Ref . |27| , and it is considerably smaller than the results obtained 
by applying the simple vector meson dominance. The SU(3) symmetry breaking in (rp) 1 / 2 , 
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FIG. 5: Dependences of {r 2 B+ B o B ) on m c = 4.2 ~ 5.2 GeV. 



TABLE IV: Mean square radius (rp) (fm 2 ) of this work and the other theoretical calculations. 



D+ 
D° 

Df 
B+ 
B° 
B° 

8 

B+ 



This work 



0.165 



-0.010 
+0.011 



-0.261 
0.0902 
0.273 
-0.134 



-0.018 
-0.019 

-0.011 
+0.014 

-0.043 
0.059 

f0.022 
-0.029 



-0 0676+ ' 0141 

u.utwo_ 0189 
0.0277 ~ 0.0451 



Lattice [27] 



0.334 ± 0.003 



VMD 



0.393 a 



"This value is obtained by (?" 2 )vmd = 6/M 2 . 
b This value is obtained by f Bc = 360 - 517 MeV. 



which mainly come from the mass difference m s — m q = 180 ± 50 MeV, are obtained as 



\ (rl + ) 



0.740 



-0.041 
+0.050' 



\ (r 2 B o) 



q 7 i 1-0.049 
u - ' ±1 +0.058' 



The radius ratio of Bt and B + are also obtained as 



0.407" 



-0.037 
-0.038 



0.319 



hO.029 
-0.030' 
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which corresponds to the range fu c = 360 ~ 517 MeV. 

V. CONCLUSIONS 

In this study, we discussed the ratios of decay constants and mean square radii for pseu- 
doscalar heavy mesons. By considering the hyperfine interaction inside the meson, we found 
that the ratio of light quark masses m s j m q was the important factor for determining the ratio 
of the decay constants. First, in the nonrelativistic approximation, we obtained the relation 
JbJ Jb = (l-05±0.02)/£i s // o which did not use any parameters. These two ratios were indi- 
vidually evaluated by including the "best-fit" light quark masses m s /m q = 483/310 = 1.558 
and the values f D JfD = 1-226 ± 0.002 and f B J Ib = 1-24 ± 0.02 were obtained. Sec- 
ond, in the light-front framework, we utilized the mass difference of light quark masses 
m s — m q — 180 ± 50 MeV and the fittings of the decay constants for light mesons to com- 
pare the mean square radii of K +, ° mesons in the power-law and Gaussian momentum 
distribution amplitudes. The latter was consistent with the data and it extracted the light 
quark masses ratio m s /m q = 1.72 ± 0.20. This mass ratio led to /d s //d = 1-29 ± 0.07 
and JbJ Jb = 1-32 ± 0.08. The former was in agreement with the experimental data and 
the latter was almost larger than all other theoretical calculations. Both these ratios were 
not only insensitive to the heavy quark masses m c ^, but also insensitive to the decay con- 
stants fo,B- Similar to the above, the ratio Jb^I Jb was mainly determined by the mass 
ratio m c /m q and the mass splitting AMb c b*- The dependences of f Bc / ' Jb on AMb c b* with 
the varied charm quark masses have been shown. We found that f Bc / 1 Jb = 1-76 ~ 2.53 
corresponded to AM Bc b* = 16 ~ 64 MeV. In addition, the mean square radii of heavy 
meson were estimated. We found the mean square radii (r 2 D+ D0 D ) and (r B + B0 Ba ) were 
insensitive to the heavy quark masses m c and m^, respectively, which was consistent with 
the behavior when the heavy quark limit was taken. Our (r B+ ) was slightly lower but 
comparable to that of lattice calculation, and was considerably smaller than that of vector 
meson dominance (VMD). The light quark mass ratio and the range of f Bc given above 
also led the radius ratios yj (r 2 Dt ) / (r 2 D+ ) = 0.740;°;^, ^/(r|o)/(r| ) = 0.7117.°;°^, and 
^/(r 2 Bt )/(r B+ ) = 0.407±HI ~ 0.319^;°!o, respectively. 
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